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We develop a theoretical approach for the dynamics of Rydberg excitations in ultracold gases, 
with a realistically large number of atoms. We rely on the reduction of the single-atom Bloch 
equations to rate equations, which is possible under various experimentally relevant conditions. 
Here, we explicitly refer to a two-step excitation-scheme. We discuss the conditions under which our 
approach is valid by comparing the results with the solution of the exact quantum master equation 
for two interacting atoms. Concerning the emergence of an excitation blockade in a Rydberg gas, our 
results are in qualitative agreement with experiment. Possible sources of quantitative discrepancy 
are carefully examined. Based on the two-step excitation scheme, we predict the occurrence of an 
antiblockade effect and propose possible ways to detect this excitation enhancement experimentally 
in an optical lattice as well as in the gas phase. 
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I. INTRODUCTION 

The possibility to routinely create samples of ultra- 
cold gases in the /i-Kelvin regime has opened a new av- 
enue to the investigation of interacting many-particle sys- 
tems. At such temperatures, the thermal velocities of the 
atoms are low enough that the atoms move a negligible 
distance over the duration of the experiment. Hence, 
(thermal) collisions are not relevant and it is possible to 
study quasi-static interactions between the particles. 

For densities of a dilute ultracold but non-degenerate 
gas typical for atoms in magneto-optical traps, the in- 
teraction between ground state atoms is very weak. Ry- 
dberg atoms, on the other hand, can strongly interact 
among each other, even in a dilute gas, due to their large 
polarizability which scales with the principal quantum 
number n ds. rJ . This scaling allows an accurate control 
over their interactions [ij over a huge range by varying 
n. In contrast to amorphous solids, with which ultra- 
cold Rydberg gases share some similarities, the atoms are 
practically stationary on the timescale of electronic dy- 
namics because of their low thermal kinetic energy 0, Q . 

A striking consequence of the strong Rydberg- Rydberg 
interaction is the so-called "dipole blockade", i.e., a sup- 
pression of Rydberg excitations due to an induced dipole 
coupling of the Rydberg atoms to their environment. 
This phenomenon was first considered theoretically in 
proposals to build fast quantum logic gates Q , to im- 
prove the resolution of atomic clocks l5| and to create 
single-atom and single- photon sources [6|. It was exper- 
imentally verified for second-order dipole-dipole (or van 
der Waals) coupling between the Rydberg atoms 0, [1] by 
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measuring the density of the Rydberg atoms as a function 
of increasing laser intensity, atomic density or principal 
quantum number, i.e., as a function of increasing inter- 
action strength. By applying and varying an external 
electric field the blockade effect was also demonstrated 
for a direct (i.e. first order) dipole-dipole interaction of 
the Rydberg atoms and it was shown that the suppres- 
sion of excitations is particularly pronounced at the so 
called Forster resonances Q. Furthermore, it was shown 
that the blockade effect also leads to a quenching of the 
probability distribution for the excitation to the Rydberg 
state 

The theoretical description of this laser-driven, inter- 
acting many-particle system is challenging. In 7] a mean 
field approach was used and the Bloch equations for a 
single Rydberg atom in a sphere were solved. Within 
the sphere, embedded in a constant background density 
of Rydberg atoms, no further excitations were allowed. 
With the help of a fit parameter the experimental results 
of |7|] could be reproduced. 

The system was also investigated by solving the many- 
particle Schrodinger equation numerically [13|. There, 
intelligent use was made of the fact that the blockade it- 
self reduces the number of atoms which can be excited 
which allows a substantial reduction in the number of 
states that had to be considered for the calculations. 
Yet, the number of atoms that could be simulated was 
still so small that appropriate boundary conditions had 
to be used to establish contact with the experiments. 
However, experiments using a two-step (three-level) ex- 
citation scheme could not be described since important 
effects, such as radiative decay, were not included. 

Here, we focus in particular on the two-step excitation 
scheme, used in the experiments (sl. [l0|. where the inter- 
mediate level decays radiatively. As we will show, this 
leads to a reduction of the description of the Rydberg 
excitation dynamics in a single atom to a rate equation 
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which in turn enables us to formulate a quasi-classical 
approach taking fully into account all atoms in the exci- 
tation volume and all interactions of the Rydberg atoms. 

Experimentally, a gas of atoms is prepared in a 
magneto-optical trap (MOT) with peak densities up to 
10^^ cm~'^ at temperatures of about 100/iK. Under these 
conditions the gas is far from the quantum degenerate 
regime and can be viewed as a classical ideal gas. Fur- 
thermore, the laser intensities used in [1] and [l^ are 
relatively low, so that coherent coupling of the atoms 
by the laser field, e.g., through stimulated emission and 
reabsorption of photons, is negligible. However, the in- 
teraction of the individual atoms with the laser fields has 
to be treated quantum mechanically. 

Our approach is based on the observation that, under 
the conditions of the experiments @ and [l^, the de- 
scription of the single-atom excitation dynamics can be 
reduced substantially to a single rate equation using an 
adiabatic approximation for the coherences. Despite the 
approximations made, the rate equation accurately de- 
scribes the population dynamics of the Rydberg state, 
including non-trivial effects like the Autler-Townes split- 
ting of the excitation line. This simplification in the 
description of the single-atom dynamics is the key that 
ultimately allows us to fully account for the correlated 
many-particle dynamics with a simple Monte-Carlo sam- 
pling, thereby reducing greatly the complexity of a full 
quantum treatment. 

The paper is organized as follows. In Sec. |TT]we present 
the approach, which enables us to describe the dynamics 
in an ultracold gas of interacting three-level atoms using 
a many-body rate equation. Starting from the full quan- 
tum master equation, we justify our approximations first 
on the single-atom level (Sec. Ill A[) . then for the inter- 
acting system (Sec. Ill B[) and finally describe how the 
Rydberg-Rydberg interaction is included in our descrip- 
tion (Sec. Ill Cp . For two interacting atoms, we compare 
the results of our rate equation with the solution of the 
quantum master equation (Sec. IIIip . In Sec. llVl we com- 
pare the results of the simulations for a realistic num- 
ber of atoms with experimental data and comment on 
the possibility to experimentally observe an interaction 
induced enhancement of Rydberg excitation ( "antiblock- 
ade"). Section IVl summarizes the results. 

Throughout the paper atomic units will be used unless 
stated otherwise. 



II. TWO-STEP RYDBERG EXCITATION IN AN 
ULTRACOLD GAS 

A. Dynamics of the non-interacting system 

In what follows, we will discuss a two-step cw- 
excitation scheme for the Rydberg state (see Fig. [Ij, as 
typically used in experiments. In the first step, the atom 
is excited from its ground state \g) to an intermediate 
level \m) with a transition strength given by the Rabi 
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FIG. 1: Sketch of the two-step excitation scheme for rubid- 
ium. 



frequency Q. The photon for this step is typically pro- 
vided by the MOT lasers, which are tuned on resonance 
with the transition \g) \m) during the time of Rydberg 
excitation. 

In the second step, a separate tunable laser drives the 
transition between the intermediate level and the desired 
Rydberg state |e) with Rabi frequency lu, where in this 
step we allow for a detuning A from resonant excitation. 
The decay of the intermediate level with rate F has to be 
taken into account, as its radiative lifetime is typically 
shorter than the pulse duration. On the other hand, the 
lifetime of the Rydberg state is much longer so that its 
decay can be neglected. 

The coherent dynamics of N non-interacting three- 
level atoms coupled to the two laser fields is governed 
in the interaction picture by the Hamiltonian Hq, 



N 



where 



/i« = A|e.)(e,| , 



N 
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+ ^i\e^){m^^ + \m^^{e^\) (2b) 

describe the interaction of the levels of atom i with the 
laser beams. 

The time evolution of the system including the decay of 
the intermediate level is then given by a quantum master 
equation for the iV-particle density matrix p^^-* = p, 



dt 



p^-t[Ho,p]+C\ 



(3) 



where the spontaneous decay of level |m) is included via 
the Lindblad operator C. In general, the rate of spon- 
taneous decay of an atom is influenced by the presence 
of other atoms through a coupling mediated by the ra- 
diation field, which can account for collective effects like 
superradiance. The strength of this coupling is deter- 
mined by the dimensionless quantity Xij = 2n\ri — rjl/A, 
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which measures the atom-atom distance in units of the 
wavelength A of the \g) — > \m) transition. For Xij <^ 1 
the spontaneous decay of an atom is strongly affected by 
its neighbors, while for Xij 3> 1 the atoms radiate inde- 
pendently. In typical experiments with ultracold gases, 
the mean atomic distance between atoms is a ^ 5/im. 
For the 5s —>■ bp transition of Rb this corresponds to 
Xij ~ 40. Hence, the collective decay is negligible and 
the Lindblad operator can be written as a sum of single- 
atom operators, 



N 

E 



with 



Li ^ \gi){mi\ and if = \mi){g, 



(4) 



(5) 



Hence, the dynamics of the atoms is completely decou- 
pled and the A'^-atom density matrix factorizes, p — 
® . ■ ■ ® p^^\ The time evolution of a non-interacting 
gas of three-level atoms is therefore completely deter- 
mined by the master equation for the single-atom density 
matrix = p, i.e., the optical Bloch equations (OBE) 
for a three-level atom. 
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As usual, the level populations are described by the di- 
agonal elements of the density matrix, whereas the off- 
diagonal elements, i.e., the coherences, contain the in- 
formation about the transition amplitudes between the 
levels. Conservation of probability leads to the sum rule 



Paa = 1 



for the populations so that 8 independent variables re- 
main to be solved for. 

This single-atom description is too complex to serve as 
the basis of a tractable description for the many-particle 



system. Fortunately, under the set of relevant experi- 
mental parameters Eq. ^ simplifies substantially. In the 
experiments [1, the upper transition is much more 
weakly driven than the lower one (w <^ ^l) due to the 
different transition dipolc moment of the respective exci- 
tation. This defines two well separated time scales, such 
that the Rydberg transition |m) — > |e) is slow compared 
to the pump transition \g) |m). Thus, the time evolu- 
tion of the system is governed by the slow Rydberg tran- 
sition in the sense that the coherences of the fast pump 
transition will adiabatically follow the slow dynamics of 
the Rydberg transition. 

Furthermore, the decay rate of the intermediate level 
is much larger than the Rabi frequency of the upper tran- 
sition (F 3> uj) implying that the populations will evolve 
only slightly over a time F"^. Hence, dephasing of the 
atomic transition dipole moments, i.e., damping of the os- 
cillations of coherences, is fast compared to the dynamics 
of the Rydberg population. 

Under these conditions, the coherences can be ex- 
pressed as a function of the populations at each instant of 
time, i.e., their dynamics can be eliminated adiabatically 
[III by setting 
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for a ^ l3 . 
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Solving the algebraic equations arising from ^ and ([S]) 
for the populations, making use of ([7]) and inserting into 
the differential equation for pmm and poc one arrives at 
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where the coefficients qk = 9^(0, a;,F, A) are some func- 
tions of the parameters describing the excitation dynam- 
ics of the three-level system. 

To simplify further, we note that within the adiabatic 
approximation ([5]) the dynamics of the population differ- 
ence Pmm — Pgg can be neglected for times t > l/2r. This 
can be verified by direct integration of pmm — Pgg from 
the OBE, which shows that the dynamics of the popula- 
tion difference is proportional to 1 — exp(— 2Ff) and thus 
reaches its saturation limit at a timescale on the order of 
1/2F. Using the sum rule ([7]) this leads to the relation 
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Pec = , 



(10) 



which can be used to eliminate the population of the 
intermediate level occurring in Finally, one arrives 
at a single differential equation for poo 



7T 



-7T : 



(7) which can readily be solved to give 



Peeit) 



1 — exp 



7T 



(11) 



(12) 



where p"^ = p^{U,uj,T,A) denotes the steady-state oc- 
cupation of level |e) and 7^ = 7|(f2,w,F, A) is the rate 
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FIG. 2: Population of the Rydberg level for the three-level 
system of Fig. [1] according to the RE Ea. (|14|l (solid lines) 
and OBE Eq.([6} (dashed lines) for different pulse lengths: 
0.3 fis (lowest pair of curves), 1.0 /is (middle pair) and 2.0 /is. 
The parameters in MHz are: (a) {n,ui,T) — (4,0.2,6) in (a) 
and (22.1,0.8,6) in (b). 



for populating the Rydberg level for short times. The ex- 
pressions for 7| and p'^ are given in the Appendix. Here 
we note only that in the limit f2 ^ F ^ w they reduce 
to 



7T 



2(1 -4A2/1]2)2 



Pc 



(13) 



which shows that the resonant excitation rate is propor- 
tional to (aj/f2)2. 

Introducing an effective ground state = 1 — pcc, 
one can write pT|) in the form of a rate equation (RE) 
for an effective two-level atom 



Pcc{t) =7T Pgf-7i Pec, 



with de-excitation rate 



7i = 7T 



(14) 



(15) 



A comparison of the solutions of the OBE ^ and the 
RE p4|) for the Rydberg populations as function of the 
detuning A is shown in Fig. [2] for different pulse lengths. 
The parameters correspond to those of the experiments 
[To| and Q . The agreement of the solutions is generally 
good and becomes even better for longer pulses. For the 
parameters of experiment [l3] the convergence of the RE 
solution to that of the OBE in the region around A = is 
slower as a function of pulse length. This is due to < F 



which indicates that it is not fully justified to neglect the 
nonlinear short-time population dynamics. 

The RE reproduces the Autler-Townes splitting of the 
intermediate level |m) manifest in a splitting of the Ry- 
dberg line, proportional to Q for short times. The split- 
ting is transient, as the steady state with its single central 
peak is approached for long times when the Rydberg pop- 
ulation reaches the saturation limit. A detailed analysis 
of the peak structure of the Rydberg populations in this 
system, especially the occurrence of the Autler-Townes 
splitting and its impact on the excitation blockade, has 
been given in [l5| . 

For future reference, we will cast the single-atom RE 
(fT4|l into a form which will be used for the simulation 
of the interacting many-particle system. To this end, 
we denote the state of the atom by a, where cr = 1 if 
the atom is in the Rydberg state, and a — otherwise. 
Furthermore, we define the rate of change for the state 



7(A,a) EE (1-a) 7t(A)+a7x(A) 



(16) 



which describes excitation of the atom if it is in the (ef- 
fective) ground state {a — 0) and de-excitation if it is 
in the excited state {a = 1). Using these definitions, we 
can combine (jl4p which determines pcc{t) and the corre- 
sponding equation for Pgg(i) in the form of an evolution 
equation for the single-atom state distribution function 



dt 



Y^T{a,o')p{a') 



(17) 



with p(0) 
matrix 



1 — Pee, = Pee and the transition rate 



T(a, a') = -7(A, a) 5„^,. + 7(A, 1 - a) 



(18) 



The first term of p8)) describes the transition a ^ 1 — a, 
through which the system can leave the state cr, while the 
opposite process (1 — a — > tr), which brings the system 
into the state a, is described by the second term. 

Proceeding to the case of N non- interacting atoms, 
we define the many-particle state cr as the configu- 
ration containing all single-atom states Ui, i.e., cr = 
{(7i, (Ti, ... (Tat} and er; as the many-body configu- 
ration which is identical to cr except for the state of 
atom i, i.e., er^ = {cti, . . . , 1 — Ci, . . . , ctat}. If we fi- 
nally use the notation 7(A, cr) = J2i 7(^, '^i) and Sa-^a-' = 
^rTi.(7[ • ■ ■•^cTN.o-^, the matrix of the transition rates gen- 
eralizes to 



T(^,^') = -7(A,<t)5„ 



(19) 



and the evolution equation for the many-body state dis- 
tribution function P{cr) can be written in a closed form 
as 



(20) 
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For non-interacting particles the rate 7 depends (besides 
on the laser detuning) only on the state of particle i, i.e., 
on Ui. However, this is no longer true in the interacting 
case and 7 will depend on the entire many-body config- 
uration. 



B. Correlated many-particle dynamics 

In order to study the correlated dynamics of the inter- 
acting many-particle system, we have to add the Hamil- 
tonian describing the Rydberg-Rydberg interaction 



H 



RR 



- E 



Uij I Ci , (^j) {^i: ^j\ 



(21) 



to Hq (c.f. eq. (HI)), where Uij is the interaction energy 
of a pair of Rydberg atoms at a distance = 
The quantum master equation ^ then reads 



-p = -I [Ho + Hrr, p]+C [p] , (22) 



with the Lindblad operator given in (j4]). 

To see which terms of the master equation are affected 
by the inclusion of the Rydberg-Rydberg interaction we 
consider the commutator [H^ + Hfm, p] in the many- 
body basis ja) = \ai, . . . un) — ■ ■ ■ law), where \ai) 
denotes the state of atom i. 



([i?A+i?«i?,,p])„;3 = E 



(a+ ^ ^5„^_^)5„^^^_(a 



Pa/3 : 



(23) 



and rewrite it (using the conservation of probabilities for and the symmetry of the Rydberg-Rydberg interaction, 
each atom, i.e., 1 = (J^^^^g -I- (5^^^,„ -I- 8^,^^^ = (J^^^^g -I- J^^^^e, = Uj-) as 



{[H,, + HRn,p])^p = Y^{5, 



5ai,g %,e) 



Pais 
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(24) 



In the first term of ([24| the Rydberg-Rydberg inter- 
action shows up as an additional (local) detuning of an 
atom at r^, whenever the atom at Vj is in the Rydberg 
state (i.e., if aj = (3j = e). In particular, no additional 
coherences are generated by the Rydberg-Rydberg in- 
teraction and, therefore, this term does not change the 
structure of the master equation as compared to the non- 
interacting case. 

The second term describes direct transitions between 



states where atoms i and j are not in the Rydberg state 
and the state where the atoms form a Rydberg pair. 
These transitions require the simultaneous absorption or 
emission of at least two photons and are thus higher order 
processes. The dynamics of these multi-photon processes 
is very slow compared to all other transitions in the sys- 
tem, therefore it can be neglected (see also the discussion 
in sections IlIII and [IV Ap . i.e., the commutator ([M]) can 
be approximated as 



J 



{[Ha+ Hrr,p])^P « ^ (5a.,e5ft,g - ^a,,g<5ft,e) 



Pa/3 • 



(25) 



Thus, within this approximation, we recover the simple picture which is commonly used for the explanation of 
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the dipole blockade effect, namely that a highly excited 
atom shifts the Rydberg levels of nearby atoms out of 
resonance with the excitation laser. 

By neglecting multi- photon transitions, the structure 
of the master equation is not changed compared to the 
non-interacting system and we can perform the adia- 
batic approximation discussed above. Identifying finally 
5aj,e ^i3j.e with (Tj , it is Straightforward to generalize Eq. 
(fT^ to the interacting case, 

T{(T, a') =. -7(A, (t) S^^„, + 7(A., 1 - a,) b^^^^. , 

(26) 

where now 7(A,cr) = 7(A, (5^, Ci) and all atoms are 
coupled by the energetic shift caused by the Rydberg- 
Rydberg interaction 



Ai = A + 5i = A - 



(27) 



so that in the interacting case the rate for a state change 
7(A,(5i,(Ti) for the atom i depends on the entire many- 
body configuration through the local detuning bi. 

The above approximations simplify the description of 
the correlated many-particle dynamics to a high degree, 
since a particular many-particle configuration cr is di- 
rectly coupled to " only" iV configurations cr' by the tran- 
sition rate matrix T(cr, <t'), which has to be compared to 
the available number of 2^ many-particle states. To ex- 
plicitly show this simplification, we insert (|26p into the 
evolution equation (|20p of the state distribution function, 
perform the sum over cr' and finally arrive at 

i i 

(28) 

Knowing Uij, Eq. (|28p can be solved with standard 
Monte-Carlo sampling techniques, allowing us to treat 
systems up to several 10^ atoms. 

We emphasize that the description presented above is 
not restricted to the three-level scheme considered in this 
work. It can, e.g., also be applied for a direct excitation 
of the Rydberg state from the ground state (two-level 
scheme) provided that the atomic coherences are damped 
out fast enough to not significantly affect the population 
dynamics of the Rydberg state (e.g., if the bandwidth of 
the excitation laser is larger than the Rabi frequency of 
the transition). For a single-step excitation scheme the 
(de)-excitation rates are given by 

2rrj2 

where F is the measured width of the excitation line. 



C. Determination of the Rydberg- Rydberg 
interaction 

An accurate determination of the interaction potential 
Uij is challenging due to the mixing of a large number of 



electronically excited molecular potential curves. Results 
from a perturbative treatment exist for the — > oo 
asymptote of the alkali-metal atoms [IB] and for the level 
shifts of Rb as well as calculations for Cs based on 
the diagonalization of the interaction Hamiltonian of two 
highly excited atoms using a large number (^ 5000) of 
pair states as basis [l^ . In the latter spirit, a simple 
picture was formulated in [l9] for Rb that allows for an 
intuitive understanding of the basic dependence of Uij 
on and on the principal quantum number n of the 
Rydberg state. 

Following [l^, a pair of Rydberg atoms in states \a) 
and 1 6) at distance experiences a shift Uij of its elec- 
tronic energy due to an induced dipole coupling Vij — 
P^aa' P^bb' I r^j to an energetically close pair of states \a') 
and 1 6'). The shift is given by the eigenvalues 



8q± 



(29) 



of the two-state Hamiltonian matrix 



where is the asymptotic {rij oo) difference between 
the energies of the two pairs. 

For a pair |ns, ns) of two atoms in the ns state, the 
relevant dipole coupling is to the energetically close pair 
\{n — l)p3/2, For an arbitrary but fixed quantum 

number tiq we may define ^^(no) = tJ'nQs(no-i)pl^nosnQp- 
The interaction strength for other Rydberg levels n then 
follows from the scaling [l[ 



M (") = M ("o) — 



6o{n) = <5o(no) ^ , 



(30a) 
(30b) 



where n* — n — rj includes the appropriate quantum de- 
fect T] (for the ns states of Rb 77 = 3.13). For rij 00 
one recovers the familiar van der Waals r~^-dependence 
and the dominant n^^ scaling for the pair interaction 
Uij. For Rb we will use in the following the values 
/i2(no) = 843800 a.u. and 6o{no) = -0.0378 a.u. for 
no = 48 from M- 



III. AN ACCURATE TREATMENT OF TWO 
INTERACTING ATOMS 

As a test for our rate equation approach in the case 
of interacting atoms, we have numerically solved the full 
quantum master equation ([^ and the rate equation (^5]) 
for two interacting atoms separated by an interatomic 
distance r. The quantity directly accessible in the exper- 
iments is the fraction of excited atoms /g. It is shown in 
Fig. [3^ and Fig.[3]D as a function of the principal quantum 
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FIG. 3: Comparison of the solutions of the master equation 
(|22|) (dashed hues) and the rate equation (I28p (solid lines) for 
two interacting atoms at distance r = 5/xm. Upper graphs 
(a,b) show the fraction of excited atoms /e, lower graphs (c,d) 
the probability pee;ee that both atoms are in the Rydberg 
state as function of the principal quantum number n for a 
pulse length r = 2 fis. The parameters of (a,c) and (b,d) are 
those of Figs. [2^ and[2]3, respectively. 



number n for excitation parameters used in the experi- 
ments and respectively. The overall agreement 
between the exact result and our approximation is very 
good and the discrepancy of only a few percent between 
the solutions is comparable to that of the single-atom 
calculations (c.f. Fig. [21 note the different scaling of the 
ordinate) and practically independent of the interaction 
strength. This indicates that most of the deviation is a 
consequence of the approximations already introduced at 
the single-atom level. 

For both parameter sets we see a suppression in /e for 
large n, i.e., an excitation blockade. Additionally, in the 
case where the single-atom excitation spectrum shows a 
double-peak structure [Fig.Oa], there is an excitation en- 
hancement for a certain n. Its actual value depends on 
the separation r of the atoms, so that in a gas this "an- 
tiblockade" will be smeared out due to the wide distri- 
bution of mutual atomic distances. However, for atoms 
regularly arranged in space, i.e., on a lattice where the 
interatomic distances are fixed, the antiblockade should 
be clearly visible To verify that the observed (anti- 
)blockade in /e is really a suppression (enhancement) of 
Rydberg pairs we have plotted the probability Peeiee that 
both atoms are in the Rydberg state. Indeed, we observe 
a complete suppression of the pair state in the blockade 
regime (Fig. [3J:) and the antiblockade peak (Fig. [3Ji) as 
well as a good agreement between the solutions of the 
master and the rate equation in both cases. 

Neglecting two-photon transitions (the second term in 
Eq. (|24p ) is the central approximation which we make in 
the description of the dynamics of the interacting sys- 
tem. In fact, these processes can be dominant, if the 
two-photon detuning vanishes far away from resonance, 



i.e., if A2ph = 2A + U{r) = for |A| > 0, \U{r)\ > 0. 
This is clearly seen in Fig. 3^, where pee;ee is shown as 
a function of the laser detuning A for two atoms sepa- 
rated by r = 5 /im. The solution of the master equation 
exhibits a triple peak structure with the central peak lo- 
cated at A = ~U{r)/2 (c.f. Ea. ([M)) ). which is not present 
in the solution of the rate equation. However, the prob- 
ability for this two-photon transition is too small to be 
visible in the signal of the total probability /e that the 
atoms are in the Rydberg state (see inset). 

Increasing the interatomic distance to r = 7/xm, i.e., 
decreasing the interaction strength, we expect that the 
blockade mechanism becomes ineffective and the contri- 
bution of Rydberg pairs to /e becomes relevant. This 
is indeed reflected in the fact that the peak of Pee-ee hi 
Fig. [4Jd is orders of magnitude higher than in Fig. |4^. 
Here, however, the atoms are successively excited to the 
Rydberg state by two single-photon transitions. Hence, 
the peak in pee-ee is correctly reproduced by the rate 
equation. 



IV. RYDBERG EXCITATION IN LARGE 
ENSEMBLES AND COMPARISON WITH THE 
EXPERIMENT 

A. Dipole blockade 

1. The density of Rydberg atoms 

We have calculated the density of Rydberg atoms as 
a function of the peak density of a Rb gas in a MOT 
according to Eg. ipS)) for excitations to the 625* and 825* 
state via the two-step excitation scheme as measured in 
8]. 

More specifically, we have determined the Rabi fre- 
quency of the first excitation step by using the data 
for the 5Si/2{F = 2) — » 5P3/2(-F = 3) trapping transi- 
tion of s'^Rb [23| and by taking the intensity of the MOT 
lasers from the experiment [21| . The measurment of fl as 
a function of the intensity of the MOT lasers using the 
Autler-Townes splitting of a Rydberg line f2l| is in very 
good agreement with our result. 

To obtain the coupling strength uj of the Rydberg tran- 
sition we have fitted it to the low-intensity measurements 
in [8] using our rate equation and scaled the result to high 
intensities and/or excitations to different principal quan- 
tum numbers. 

Fig. O shows the results of our calculations and the 
experiment. Although we see a qualitative agreement 
we predict Rydberg densities about twice as large as the 
measured ones. As the curves for both principal quantum 
numbers exhibit the same deviation from the measured 
data, it is tempting to scale our results to the experimen- 
tal points using a common factor. Note, however, that 
without other influences in the experiment, there is no 
free parameter in our description that would justify such 
a scaling. 
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FIG. 5: Density of Rydberg atoms as a function of the peak 
density in the MOT for a pulse length of r = 20 fis for the 82S 
(black) and the 62S state (gray) of Rb. Circles: experimen- 
tal data taken from Q|. Lines: Calculations using different 
models for the pair interactions potential: two-state model 
of ref. [31 (solid) and pure van der Waals interaction from 
perturbative treatment [Ig] (dashed) . 



cients C(i{n) for the latter are calculated in second-order 
perturbation theory for r — » oo and have been taken 
from [Tg] . The interaction strength for the nS states cal- 
culated in this "way is considerably larger than the one 
from the t"wo-state model (e.g., for the 825" state the dif- 
ference in C/(r) at r = 10 fj,m is roughly a factor of 2.5 and 
increases "with decreasing r). Yet, the final results for the 
Rydberg population differ only slightly (see Fig. [5|) . We 
conclude that pc is relatively robust against changes in 
the interaction strength. This is due to the fact that the 
measurement of the Rydberg density as a function of the 
ground state density does not probe the exact shape of 
the interaction potential but rather the critical distance 
Tc at "which the energetic shift caused by the interaction 
becomes larger than half the -width of the spectral line 
(« 20 MHz). For U{r) determined in perturbation theory 
and estimated by the t"wo-state approximation Tc ~ 8 /im 
and Tc = 7 fj,m, respectively, for the 82S state, so that 
significant differences emerge only for large densities. 



The influence of ions 



In the follo"wing we estimate the quantitative influence 
"which several effects could have on the results presented. 



2. The influence of different Rydberg -Rydberg interactions 

The "exact" Rydberg-Rydberg interaction may differ 
from the one we have used in our description. To assess 
the impact of such a difference, we have performed our 
calculations with the simple two-state model discussed 
above (solid lines in Fig. [5|) and assuming a "pure" van 
der Waals interaction, —C^/r^, between the Rydberg 
atoms (dashed lines in Fig. [5]). The interaction coeffi- 



Another effect, so far not accounted for, is the pres- 
ence of ions. The excitation pulse length used in Q was 
20 fj,s. For pulse durations that long, it was shown that 
a significant amount of Rydberg atoms can undergo ion- 
izing collisions even for a repulsive Rydberg-Rydberg in- 
teraction [I^, [22|. The presence of ions in the system 
influences the excitation dynamics due to the polarizing 
effect of the electric field of the ions on the highly sus- 
ceptible Rydberg atoms. The Rydberg- ion interaction 
(oc r^*), therefore, leads to an additional energetic shift 
of the Rydberg levels and, thus, can lead to an enhanced 
excitation suppression. 

To see if the presence of ions can account for the dif- 
ference between our results and the measured data, we 
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FIG. 6: Estimated average number Np of n = 82 Rydberg 
pairs excited by multi-photon transitions as a function of the 
laser detuning A after r = 20/is for a ground state peak den- 
sity po = 10^° cm"^. 

have performed calculations in which we have replaced 
up to 20% of the Rydberg atoms by ions. The change 
in the results compared to the situation without ions is 
comparable to that of stronger Rydberg-Rydberg inter- 
action discussed above. Therefore, ions can be ruled out 
as a source for the discrepancy between our and the ex- 
perimental results. 

4-. The influence of multiphoton transitions 

The excitation line profiles presented in Q showed an 
enormous broadening for measurements at high densities. 
In contrast, the line profiles that we have calculated with 
the present approach are much narrower, in accordance 
with the simulations reported in ref . ^IS] . 

The strong line broadening in the experiment could be 
due to non-resonant effects, such as multiphoton transi- 
tions, not included in our rate description (see discus- 
sion in section IIII[) . To estimate their possible influ- 
ence, we have to determine first the number of Ryd- 
berg pairs which could be excited by these transitions. 
To this end, we have determined the number of (ground 
state) atoms np{r)Ar which form a pair with a distance 
between r and r -I- Ar in the excitation volume, from 
the pair density np{r). Furthermore, we have calculated 
the probability Pee;ee for a pair of atoms to be in the 
Rydberg state after r = 20/is by solving the quantum 
master equation (Wj^^) and the rate equation [wj^^) 
for two atoms as a function of the laser detuning A 
and interatomic distance r (c.f. Fig. 0]). The difference 
Wp{r, A) — Wj^^{r, A) — Wj^^{r, A) should give a rough 
estimate for the probability of a Rydberg pair being ex- 
cited by a multi-photon transition. The average number 
of such pairs as a function of A can then be estimated 
by Np{A) = J2^Wpiri,A)np{ri)Ar. 

Fig. [5] shows that for a sample with ground state peak 
density po = 10^° cm~'^ our estimate yields a negligible 



number of Rydberg pairs excited by multi-photon tran- 
sitions after 20/is. Although these estimates are rather 
crude, the result shows that multi-photon effects are too 
small to explain the broadening of the excitation line pro- 
file in the experiment @ . 

In summary, the unexplained line broadening and the 
difference between experiment and theory in the Rydberg 
populations make it likely that some additional, presently 
not known process, has contributed significantly to the 
results obtained in 

B. Antiblockade 

1. Lattice configurations 

The discussion in Sec. Illll has shown that the structure 
of the single-atom excitation line strongly influences the 
excitation dynamics in the interacting system. Even on 
resonance, the Rydberg-Rydberg interaction can cause 
an excitation enhancement, if the spectral line exhibits a 
double peak structure. This antiblockade occurs when- 
ever the interaction-induced energetic shift A^ for an 
atom at position matches the detuning A^ax at which 
the single-atom excitation probability has its maximum 
value. 

In the gas phase, where the mutual atomic distances 
are broadly distributed, the antiblockade can hardly be 
observed by measuring the fraction of excited atoms /e, 
as the condition A^ = Amax is only met by relatively 
few atoms ^15 ]. In contrast, if the atoms are regularly 
arranged in space, e.g., with the help of an optical lattice 
produced by CO2 lasers [2^, one should clearly observe 
peaks in /e for certain n (see Fig. [7^). The peak posi- 
tions can easily be determined by analyzing the geometry 
of the underlying lattice. Moreover, the effect is quite ro- 
bust against lattice defects (unoccupied lattice sites) and 
should therefore be experimentally realizable. A more 
detailed discussion can be found in [15,] . 

The underlying lattice structure allows for a statisti- 
cal interpretation of the antiblockade as "clustering" of 
Rydberg atoms. Using the terminology of percolation 
theory, we define a cluster of size s as group of s nearest 
neighbor sites occupied by Rydberg atoms. For negligi- 
ble Rydberg-Rydberg interaction the excitation of atoms 
on a lattice is analogous to the situation encountered in 
classical (site-)percolation theory. This is seen in Fig.lJb, 
where a histogram of the average number ris of s-clusters 
per lattice site as function of the cluster size (normalized 
to the number of 1-clusters, i.e., isolated Rydberg atoms) 
is shown for atoms excited to the state n = 40. The 
shaded area represents the prediction of percolation the- 
ory 2^ for the same number of isolated Rydberg atoms 
per site and shows good agreement with the "measured" 
data. In the antiblockade regime (n = 65, Fig.[7l:) we ob- 
serve a broadening of the cluster size distribution and a 
significant enhancement of larger Rydberg clusters, while 
in the blockade regime {n = 68, Fig. (TJi) a quenching of 
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FIG. 7: (a) Fraction of excited atoms for atoms on a simple cubic lattice with 20% unoccupied sites as function of the principal 
quantum number n. The lattice constant is a = 5/im, all other parameters are those of Fig.[2]D. (b-d) Corresponding number of 
"Rydberg clusters" per lattice site Us normalized to the number of 1-clusters (i.e. isolated Rydberg atoms) ni as a function of 
the cluster size s for principal quantum number n = 40 (b), n = 65 (c) and n = 68 (d). The shaded areas represent predictions 
from percolation theory 2^ for ^ system with the same number of isolated Rydberg atoms (1-clusters) per lattice site. 




Mandel's Q-parameter 



FIG. 8: Comparison of the Q-parameter in the blockade 
(squares) and antiblockade (circles) configuration as a func- 
tion of the principal quantum number n for a sample with a 
homogeneous atomic density po = 8 x 10^ cm""^ and for an 
excitation pulse length r = 2 ^s. Q was determined by 10^ 
successive measurements of (Nryd) and (N^yd)- The Rabi fre- 
quencies {n, Lu) are: (4.0, 0.24) MHz (squares) and (22.1, 0.8) 
MHz (circles). 



the distribution and an enhancement of the probability 
to excite isolated Rydberg atoms is evident. 



2. Random gases 

Based on the solution of a many-body rate equation 
using Monte Carlo sampling, the present approach is par- 
ticularly well suited to determine statistical properties of 
interacting Rydberg gases. 

In [l3| the distribution of the number of Rydberg 
atoms was measured as function of the interaction 
strength. The distributions obtained were quantified by 



(31) 



where is the number of Rydberg atoms and (...) de- 
notes the average over the probability distribution. The 
Q-parameter measures the deviation of a probability dis- 
tribution from a Poissonian, for which it is zero, whereas 
for a super-(sub-)Poissonian it is positive (negative). The 
experiment showed a quenching of the Rydberg number 
distribution, i.e., a decrease of Q, for increas ing inter- 
action strength as theoretically confirmed (Til . [l3 | . The 
differences between the theoretical calculations {Q < 0, 
for all n) and the measured values {Q > 0) can be at- 
tributed to shot-to-shot fluctuations of the number of 
ground state atoms in the experiment (25| . 

The excitation parameters in [l3| were in the blockade 
regime, where the single-atom excitation line exhibits a 
single peak at A = 0. Therefore, there is a volume ("cor- 
relation hole") around each Rydberg atom, where the ex- 
citation of additional atoms is strongly suppressed. On 
the other hand, in the parameter regime of the antiblock- 
ade, where the excitation line shows a double peak struc- 
ture, there is in addition a shell around each Rydberg 
atom, in which additional excitations are strongly en- 
hanced. Thus, the statistics of the Rydberg excitations 
should depend on the structure of the single-atom excita- 
tion line and the antiblockade can be detected indirectly 
even in the gas phase by measuring the atom counting 
statistics. 

Figure [5] shows the calculated Q-parameter as a func- 
tion of the principal quantum number n for the block- 
ade and antiblockade regime. In the blockade configu- 
ration (squares) one observes a monotonic decrease of 
Q with n in accordance with the measurements in p^ . 
In the antiblockade regime (circles), however, Q is non- 
monotonic, i.e., the distribution is slightly broadened. 



11 



and the quenching starts at much higher n. Although 
the broadening of the distribution may be difhcult to 
observe experimentally, the difference in the functional 
form of Q{n) provides a clear experimental signature in 
a mesoscopic region of the MOT, where the atomic den- 
sity is approximately homogeneous. 



CONCLUSIONS 



scheme has lead to the prediction of an antiblockade ef- 
fect due to an Autler-Townes splitting of the intermedi- 
ate level probed by the Rydberg transition in the appro- 
priate parameter regime. This antiblockade should be 
directly observable for a lattice gas, realized, e.g., with 
an optical lattice. As we have demonstrated, it could 
also be observed indirectly in the gas phase through the 
atom counting statistics which differs qualitatively from 
its counterpart in the blockade regime. 



We have developed a simple approach, which allows 
one to describe the dynamics in ultracold gases, in which 
Rydberg atoms are excited via a resonant two-step tran- 
sition. Starting from a quantum master equation, which 
incorporates the full dynamics of an interacting gas of 
three-level atoms, we have derived a many-body rate 
equation. It covers the correlated dynamics of the sys- 
tem, yet, it can easily be solved by Monte Carlo sampling 
for a realistically large number of atoms. 

Our approach, valid under well defined conditions typ- 
ical for experiments, is based upon two approximations: 
(i) an adiadabtic approximation on the single-atom level 
to eliminate the atomic coherences and (ii) the negli- 
gence of multi-photon transitions in the interacting sys- 
tem. Solving the problem of two interacting atoms ex- 
actly with a quantum master equation we could show 
that the approximate solution based on the rate equa- 
tion is in very good agreement with the exact result. 

The present approach is capable of reproducing the 
partial excitation blockade observed in |8i| qualitatively. 

Qualitatively in accordance with our calculations re- 
garding the excitation line shape and the so called Q- 
parameter are also the experimental results of p^ . 

Finally, the careful analysis of the two-step excitation 



APPENDIX: EXPRESSIONS FOR p2l AND 7^ 

The steady state solution of the OBE for the Ryd- 
berg population is 

(A.l) 



(f72 + a>2)2 + 4A2 (r2 2f72) ■ 
The excitation rate in (|11|1 can be written as 



7T 



where 



ao + a2A2 + a4A* 



00 = 



(uj'^ - 2n'^Y • "'T-2 /o2 



2T^ {n^ + uj- 
+4lu^ (r^ + 4fi2) + 8lu^ 

04 = 32 (r^ + 2n^) . 



(A.2) 



(A.3a) 

(A.3b) 
(A.3c) 
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